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I consider the effect of the finite width of the resonant level on its thermodynamics 
and tunneling transport properties in the single electron charging regime. The finite 
width of the levels results from their derealization with formation of a narrow band 
due to their mutual overlap. Making use of 1/N technique one can show that there 
is a fixed energy level position at which the system undergoes the transition from 
the Kondo regime (above this energy) to the paramagnetic regime (below it). The 
linear bias conductance calculated as a function of the chemical potential turns out 
to be suppressed all around except for the vicinity of the transition energy where it 
has an asymmetrical resonant behavior. 
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There is much interest in the study of low temperature transport in mesoscopic systems 
where new phenomena associated with quantum coherence are invoked by the large Coulomb 
charging energy of a localized constriction. These have been mostly discussed for the case of 
tunneling through a single spin-degenerate level specified with Anderson model Q. Another 
way of manifesting a system with such a quantum correlation was recently suggested for a 
conducting granule in the Coulomb blockade condition [§], f5[. In these systems, in partic- 
ular, the quasiclassical size of the granule results in the derealization of the states and the 
appearance of a continuum band instead of the separate levels. 

In this paper I will consider the effect of the derealization of the resonant states on 
the Kondo resonance phenomenon dealing with the problem of tunneling through a double 
barrier (2D) structure (DBS) located between two conductors. This structure consists of the 
levels located at the impurity sites between the barriers. The sites are ordered into a lattice, 
and the levels of the neighbor cites have a small overlap, which leads to the appearance of 
a narrow band of width Eq. First, one could discuss this system as a vertical quantum dot 
structure where the tunneling layer is inside a constriction of radius R. This band should 
be considered continuous if its level splitting A G ps hv/R ( v is the speed of the band 
electron with energy near the bottom ) is less than the inverse time T = %/r of escape of 
the band electron into the electrodes. To make it more clear one could assume that this 
layer is infinite and the effective tunneling area through the layer is only of R 2 size. Then 
the probability for the tunneling electron to escape from this area along the layer could be 
estimated as relation between the squares of appropriate surfaces 2ttRvt and ttR 2 . This 
probability becomes negligible if vr <C R. 

The charging energy for accumulation of the charge Q inside the barriers could be eval- 
uated in a standard way || as Q 2 /2C and decreases as R 2 due to dependence of the ca- 
pacitance C(R) on the area of the contact. Since the charging energy of a non-equilibrium 
distributed charge should be more than the one in equilibrium the criterion of a single- 
electron charging of the layer could be formulated as R < Iq where e 2 /C(lc) > e&, e& is the 
bottom of the impurity band. To evaluate possible parameters of such a system one could 
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use as an example the typical Metal-Oxide-Metal junction, where closely located conductors 
with the tunneling area of lOOnm 2 form a condensator with m 100K charging energy. This 
value of the charging energy is substantially greater than the band width which, due to the 
previous restriction according to T, should be less than O.IK for 10 2 — 10 3 lattice sites. The 
dependence on R in this discussion could become irrelevant if the phase coherence length of 
the electrodes is less than R. I could generalize the consideration if I will suppose that this 
length satisfies the above two criteria instead of R. 

Tunneling through the DBS is locally described by the Hamiltonian 

£a,kCk S) a C ks,a + ^ e qd^ s dqs + ^ (Va,q,kd^ s c ks,a + ^a,q,k C ks,a^qs) (1) 
k <7 k,q,a,s 

with a restriction on the Hilbert space that only states of zero |0 > and one-electron \qs > 
occupation of the DBS band are permissible. Here Ck s , a annihilates the electron of a lead 
with spin s in the state of wave vector k = (k±,k\\) and of energy e a ,k- The components 
k±, k\\ are perpendicular and parallel to the layer, respectively. The operator d qs annihilates 
an electron in the state localized inside the layer with spin s, and wave vector q, and energy 
e q . An expression for the tunneling matrix elements V aA ^ could be generally written making 
use of the tight binding model for the DBS band as 

V a ^k = * 1/2 ^2T aj (k±)exp(iRjq - i(f> aj (k\\)) (2) 
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Here M is the effective number of the sites of the layer restricted either by the geometry or 
by the coherence length. L is the normalization in the direction perpendicular to the layer, 
and T a j is the modulus of the corresponding tunneling amplitude. The vertex function may 
be introduced as r au/ (q 1 , q 2 ) = nJ2k ~ e k)Va,qi,kYa,q2,k which is completely ample for our 
future needs. Moreover, in the low temperature limit only the transitions between the holes 
of the Fermi surface (u ~ 0) and the electrons of the bottom of the layer band are important. 
( I will suppose that the bottom is reached at e q=0 — e&.) To calculate the function in this 
limit one could exploit either the jelly model for the leads or the tight binding model if all 
structure are commensurate. Virtually, the result is the same in the above limit. Indeed, 
taking <pj a = Rjk\\ one finds 
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r a , w (gi,g 2 ) = 7T ^ fopn+Qfopqa+Q' / de k± p 1D (e k± )5(u - e k± - e fe|| )T a 2 (e fc J ~ 7iT^p a)1D (0)S, 
Q,Q' 

(3) 

where Q, Q' are the vectors of the reciprocal lattice either of the DBS or of the leads. T a j is 
taken as a constant if the position of the j DBS site corresponds to the position of the lead 
site ( and in the case of the jelly model). Otherwise, it was put zero. The second equation 
in (HD is correct if the DBS lattice is less close than the lattice of the leads. If this condition 
is not met, the transitions between nonequal (71,52 become permissible. I will not discuss it 
specially, though, the consideration could be easy generalized to account for them. 

Conservation of the q components in (|j) makes the tunneling resonant and justifies the 
application of the 1/N expansion technique [p]]. In this case it is a 1/M expansion and it 
means the summing up of the most divergent terms in M of each order of expansion in T. 
I will restrict my consideration to leading order, giving summation of all terms of (TM) n . 

Starting from the thermodynamics I will include the DBS contribution into the partition 
function in terms of the self energy of the slave-boson T, (z) []]]. To leading order it takes 
the form 

Z(1) £e -^ + / jj^-^W (4) 



where C encircles all the singularities of the integrand, (5 = 1/T is the inverse temperature, 
and Eg 1 " 1 (z) is given by 

E« W = £ f ^ 7 IM- (5) 

In the low temperature limit the Fermi-Dirac function f(u) can be approximate by the step 
function. Then for 2D band, which has a constant density of states Pid — m * / (27ih 2 ) 
near the edge, is 

(z) = Ae b + P[E G ln h (e 6 - z)ln J (6) 

e b — z 

where P = 2(T/ti)p2dS ~ (T/ti)M/Eg, T = J2 a ^a, S is the square of the tunneling surface, 
Aeb — PEq{Iti{Eq/ Ep) — 1). P is a number of the states lying inside the energy interval 



of width T/tt. In the above condition P> 1. This function coincides with the one of a 
single level Kondo model for \e b — z\ ^> E G . However, for small |e& — z\ its behavior is quite 
different. In particular, the self-energy Eg (z) converges to the finite negative value Ae b at 
e b and acquires an imaginary part on the real axis above e b as ImT$(x — iO) = nP(x — e&). 

This means that unless e b becomes lower then Ae&, the lowest singularity of the integrand 
in (f|) is a simple pole z = e b — T K , which corresponds to the bound state of the slave-boson 
with energy Tk = E G x 

"' - ln(l + x) + xln(l + x~ x ) + I-, e b = e b -PE G (ln(E G /E F )-l) (7) 



PEq ' P 

In the limit e b /\Ae b \ <C 1 the solution of (^) loses its dependence on e b and shows that 
Tr ~ |Ae&| when .Bp 3> -Eg- I n the more general case e b could reach the zero, where the 
transition into a paramagnetic state of the DBS occurs. In that regime the low temperature 
thermodynamics of the DBS is ruled out by the lowest states of its band. 

The standard way of calculating || the localized charge the charge xc an d the spin 
Xs susceptibilities results, for the Kondo regime behavior in 
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dF DB s 1 
e b - 1 K , n G = — = 1 



PEgZl = PE G Z 

XC ~ T K (E G + T K ) ' Xs 8T K (E G + T K ) U 
1 1 



1 - d z Z (z)\ z=eb _ TK 1 + Pln[(E G + T K )/T K \ 
The temperature dependence is dropped in Eq. (^j). It remains unimportant while PT/T <C 
exp(Tx/T). These solutions (|7|, ||) predict the singular behavior of all thermodynamic 
characteristics at e b = 0. The function n G goes to 1 as l/ln(e b / E G ). 

In the paramagnetic regime the low temperature asymptotic for (|j) is 

Z« MT 



Z ba nd E G 



CO 

dx 



e-^ T a(e b ,T) (9) 

£ e b + Px 

7T (e b + Px(l + ln(E G /x))) 2 + (nPx) 2 

If — e b > TPln(E G /T) the parameter a is evaluated as a = 1 — T/(7re b ). From Eq. one 

concludes that n G = 1 and Xs diverges as 1/T. 



a = 1 -p ( 
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Discussing the transport through this structure I will restrict consideration to linear 
conductance. Then, the conductance can be expressed in terms of the imaginary part of the 
retarded Green function [0] of the DBS electrons G q (u) labelled with the quasi momentum 
q and energy ui indices. Up to T 2 corrections it is given by 

2e 2 r B r L 



0o 



■P£/ de g ImG>)Uo (10) 



h r 2 

In this technique the full Green function of electron is calculated as a convolution [|I|] of 
the slave-boson and the slave-fermion Green functions along the appropriate contour in the 
complex energy plane. For the imaginary part it is simplified to 

ImG,(u) + ) = - 1 I" ± I^)l m ^ {x + . ) + x 

Z/ZbandJoo n \x - h {X)\ 2 \x + UJ - e q - T, q (x + u)\ 2 
Here S g (x + ) is the self energy of the slave fermion which is taken at the real point x + = x+iO 
with the infinitely small imaginary shift. 

In the Kondo regime the calculations straightforwardly follow the ones of the single level 
case. To leading order in the 1/M expansion, the imaginary parts of both the slave boson and 
the slave fermion Green functions become strongly localized as Z 5(uj — T#) and S(u — e q ), 
respectively, at low temperatures by the factor exp(— T/Tk)- To obtain a non-zero ImG q (0) 
I should use the next order expression for ImT lq {x) standing in the numerator of (pTTj) . This 
is written as 

Im^ q {y + ) = Im- / w = -TZ (1 - f(y - e b + T K )) (12) 

y-x-Ey(y + -x) 



for y < e&. Finally, I obtain 



r n , s FZq 4e 2 r i ,r L pz 2 e g n _. 

ImG q {uj + ) = — — -y — — -, a = — u 13 

{u - e q + e b - T K y h T 2 T K (T K + E G ) 

Taking T% from Eq.(^) and Z from Eq.(^) one finds 

2e 2 T R T L r 
a ° h r 2 7Te b ln(E G /e b ) 1 ] 

for small Tk- This shows that the resonance in the o"o dependence on the bottom energy is 

reached at e& = Ae b , and it decreases with increasing e b more slowly than the l/e& 2 tail of 



the one electron resonance. The width on this side of the resonance is less than T in the 
factor l/ln(Ea/T). Eq. (|T3| ) is applicable only far from the resonance. Near the resonance 
it could be corrected with the summation of the infinite series of the 1/M expansion as the 
NCA does. 

In the paramagnetic phase, when e b lies below Ae b , the imaginary parts of the self 
energies give a contribution into the conductance already in the leading order. However, 
this contribution is oc T since Eq^(?/) oc (y — e b ). The zero temperature part of <jq appears if 
the imaginary part of Eo standing in the numerator of (|TTD is taken in the next order. Then 
it is equal to 

= P- / de[- 1 - — J— ] / — * (15) 



vr J-e f z + e~e b - E G z + e-e b Jo ^ + e _ e ' _ £ 
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Extracting from here the imaginary part and collecting the leading terms one could find 

n 1 aeb 7i 

for T < \ef,\ / (Pln^r) . The coefficient <p is given by 

r E p de 1 ErP PEn 

J o e b - e - Eq (e b - e) r |e 6 | & F + ^ti G 

Equation ( [r6| ) describes the decay of the resonance below e b = Ae&. It reveals that the 

square of the halfwidth is larger than T 2 by a factor P</>. 

Finally, one could conclude that the derealization of the resonant levels with formation 

of a narrow band which is located inside the barrier in between two electrodes, leads to the 

suppression of the Kondo conductance all around except for a special position of the level 

energies at e b = PEc^ln^Ec/ 1 Ep) — 1). At this point the tunneling structure contributions 

to the thermodynamic functions become singular, and the linear conductance has a resonant 

behavior. These results are easy related to the known properties of the Coulomb blockade 

transport through the granule . The above consideration describes the first steps of the 

Coulomb staircase. Due to the lack of the electron-hole symmetry in the excitation spectrum 

of the granule the position of the resonance is shifted from no = 1/2 to no = 1, and the 

form of the resonance undergoes asymmetrical deformation ||. 
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